Fluctuation-dissipation ratio in lattice-gas models with kinetic constraints 
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We investigate by Montecarlo simulation the linear response function of three dimensional structural 
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Introduction. Glassy materials are a paradigmatic ex- 
ample of systems characterized by slow relaxation and 
physical aging phenomena While thermodynamic 

observables are rather regular along the glass transforma- 
tion, dynamical properties undergo drastic changes, and 
in particular relaxation times increases up to 15 orders 
of magnitude. The microscopic origin of such a dramatic 
slowing down is the so called cage effect: a tagged par- 
ticle in a liquid does not perform a simple random walk 
but rattles in a cage formed by the surrounding particles, 
and moves away from it only if the particles building up 
the cage walls are themselves able to move. At low tem- 
perature or high density, when the degrees of freedom of 
the particles become strongly coupled, this cooperative 
process exceeds observation time and the glass appears 
like a "solid" @. 

However, structural relaxation is not completely ab- 
sent in glasses, and molecular rearrangements still oc- 
cur on longer time scales. Indeed, after a quench in the 
low temperature or high density phase, particle diffusion 
does not stop altogether, but rather becomes slower and 
slower as time goes on. This is the origin of aging effect: 
the response to an external perturbation depends on the 
"age" of the sample, i.e., on the time elapsed since its 
preparation, even at very long times. Therefore in this 
regime glassy systems never reach thermodynamic equi- 
librium and, consequently, time translational invariancc 
(TTI) and the fluctuation-dissipation theorem (FDT) are 
violated. 

The problem. Despite the remarkable progress real- 
ized in the study of out-of-equilibrium dynamics of aging 
systems [|| , the detailed nature of ergodicity breaking in 
finite-dimensional glasses is not yet fully understood. In 
mean-field mode-coupling theory the glass transition ap- 
pears as a purely dynamic effect due to an instability of 
the equation governing the correlation of density fluctu- 
ations |^]. In particular, mean-field disordered models of 
structural glasses show that the origin of this dynamical 



transition is the existence of a large number of metastable 
states which trap the system for an infinite time. On the 
other hand, the life-time of metastable states in finite- 
dimensional short range models is finite, since it is always 
possible to nucleate, by a thermally activated process, 
a droplet of the stable phase. Therefore the dynamical 
transition appear as an artifact of the mean-field approx- 
imation, and in real glasses this transition would be a 
finite-time effect, at least on time scales much smaller 
than the life-time of metastable states || |j| . 

However, since salient features of glassy behavior are 
essentially of dynamical nature, and since dynamical uni- 
versality classes are smaller than the static ones, it is 
important to establish how glassy effects depend on the 
details of dynamics. In this paper we explore the limit 
case of three dimensional lattice gas models defined only 
by short-range kinetic constraints and by a trivial equi- 
librium measure, in which the glass transition appears to 
have a purely dynamical origin unrelated to the existence 
of metastable states p|,p!o[. Previous Montecarlo simu- 
lations have shown that this model reproduces qualita- 
tively the "glassy" phenomenology, including annealing- 
rate dependence, irreversibility effects, and waiting-time 
dependent logarithmic diffusion Jl(],[ll|]. Here we investi- 
gate aging effects in the linear response function in order 
to characterize the violation of FDT 12 13f|. 

The model. Our starting point is the kinetic lattice- 
gas model introduced by Kob and Andersen in order to 
mimic the cage effect in supercooled liquids ||. The 
system consists of N particles in a cubic lattice of size 
L , with periodic boundary conditions. There can be at 
most one particle per site. Apart from this hard-core 
constraint there are no other static interactions among 
the particles. At each time step a particle and one of its 
neighbouring sites are chosen at random. The particle 
moves if the three following conditions are all met: 

1. the neighbouring site is empty; 
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2. the particle has less than 4 nearest neighbours; 

3. the particle will have less than 4 nearest neighbours 
after it has moved. 

The rule is symmetric in time, detailed balance is satis- 
fied and the allowed configurations have the same weight 
in equilibrium. With this simple definition one can pro- 
ceed to study the dynamical behavior of the model at 
equilibrium. One observes that the dynamics becomes 
slower and slower as the particle density p increases; in 
particular, the diffusion coefficient of the particles, D, 
vanishes as the density p approaches the critical value 
p c ~ 0.88, with a power law 



D(p) ~ ( Pc - pf, 



(1) 



with an exponent <f> ~ 3.f (9|. 

Since we are interested in the dynamical approach to 
the putative equilibrium state we allow the system to 
exchange particles with a reservoir characterized by a 
chemical potential p. Therefore, we alternate the ordi- 
nary diffusion sweeps with sweeps of creation/destruction 
of particles on a single layer with the following Monte- 
carlo rule: we randomly choose a site on the layer; if it is 
empty, we add a new particle; otherwise we remove the 
old particle with probability e - ^. The number of parti- 
cles is no longer fixed and the external control parameter 
is p, which plays the role of the inverse temperature. In 
this way we can prepare the system in a non equilibrium 
state by a process analogous to a quench, which is repre- 
sented by a jump in f / p from above to below 1/ p c (where 
p c is defined through the state equation p(p c ) — p c )- The 
situation becomes analogous to the canonic case in which 
one controls the temperature, and the energy endeavors 
to reach its equilibrium value. 

Linear response. Following a suggestion by J. Kur- 
chan, we compute the linear response function by apply- 
ing to our system a small random perturbation at time 
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where = ±1 independently for each coordinate a and 
particle k at position r£. The response function of the 
system at time t + t w is defined by 

'df£-ri(t + t w ) 
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a=l fc=l 



de(t w ) 



(3) 



The linear regime is probed for small enough values of the 
perturbation strength e. It is actually more convenient 
to look at the integrated response function: 



x(t + t w ,t w ) 



t+t„ 



R{t + t w ,T)e{ T )dT 



(4) 



In order to measure this quantity we prepare the system 
in a non-equilibrium state by quenching to the subcritical 
value of 1/p and letting it relax up to time t w ; at t w , we 
make a copy of the system and apply the perturbation; we 
then let evolve the two systems with the same succession 
of random numbers and measure the difference Ar£ (t + 
t w ) between the displacements that take place in the two 
systems at time t + t w . At constant field we obtain 
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t w )). 



(5) 
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We have performed a sudden quench to the subcritical 
value 1/p = 1/2.2 starting from a configuration with den- 
sity 0.75 and measured the integrated response function 
X(t + t w , t w ) vs. time t at different waiting times t w . We 
have checked that for 0.05 < e < 0.15 non-linear effects 
are absent. In the following we present the results for the 
case e = 0.1 and for a cubic lattice of size 20 3 . Figure |l| 
shows that for small values of t w TTI holds, and aging 
effects set in only for larger value of t w . To obtain the 
scaling behavior one has to consider an effective waiting 
time, t w = t w +TQ, that takes into account the relaxation 
time, ro, of the system before the quench. Indeed, figure 
|2| shows that, if one plots x(t+t w , tw) vs. t/r w , the curves 
lie roughly on top of each other and follow a simple ag- 
ing logarithmic behavior, x[t + t w ,t w ) ~ log(l + t/r w ). 
Similar results were previously found also for the mean 
square displacement [jlO 11 1. 

Fluctuation- dissipation ratio. We are now in a po- 
sition to investigate the nature of breakdown of ther- 
modynamical equilibrium. We consider the generalized 
FDT proposed by Cugliandolo and Kurchan in the frame- 
work of dynamical mean-field theory of spherical p-spin 
glass model Jl2| . In our case the relevant property is the 
off-equilibrium generalization of Einstein-Stokes relation 
among diffusion coefficient and viscosity, given by 

x{B) dB(t + t w ,t w ) = _ 2TR(t + (g) 



where 
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B(t + t w , U = ^EE (Kit + t w ) - r a k {t w )] 2 ) , (7) 
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is the mean square displacement. The fluctuation- 
dissipation ratio (FDR), X, depends on both times only 
through B and its departure from 1 is a measure of vio- 
lation of FDT. The integrated version of the generalized 
FDT allows to write (considering a constant field): 



x(t + t w ,t w ) = 



r B(t+t w ,t w ) 



2T 



X(B)dB 



(8) 



(We always write explicitly the temperature T, however 
it does not play any role.) Therefore, if FDT holds, we 
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would have a straight line with slope e/2T in the para- 
metric plot of B(t + t w ,t w ) versus x(t + t w ,t w ); while 
a deviation from this straight line indicates a failure of 
FDT. The way in which FDT breaks down plays a key 
role in aging systems since it signals the presence of a 
time-scale dependent effective temperature |p| , and this 
violation can be bounded by a quantity related to the en- 
tropy production |Q. The FDR has been measured by 
Montecarlo simulations in several models such as spin- 
glasses |L7 20 , Lennard- Jones binary mixtures p8[ , do- 
main growth models jl9[ ] , and directed polymers plj , and 
so far all the findings are consistent with the prediction 
of generalized FDT. We have carried out extensive nu- 
merical simulation to measure the FDR, X(B), in three 
dimensional lattice-gas models with kinetic constraints. 
We observe the two asymptotic regimes also found in 
the out of equilibrium dynamics of mean field models of 
structural glasses |l^] : 

• at times t smaller than t w we observe a short-time 
quasi equilibrium regime where FDT holds; 

• at larger separation times FDT breaks down with 
a constant violation factor. 

These results appear clearly in the parametric plot EH] 
of x(t + t w ,t w ) vs. B{t + t w ,t w ) where, as can be seen 
in fig. ||, the curve approaches the characteristic broken 
line. (In order to probe the asymptotic time regime we 
have chosen t w = 10 5 .) It is natural in this context to 
interpret the quasi-equilibrium regime as corresponding 
to the motion of particle inside a cage, although the size 
of cage in this model is quite small and of the order of a 
fraction of lattice spacing. After that the system crosses 
over a diffusion-cage regime where a departure from equi- 
librium is observed with a constant FDR. 

In order to see if the inverse chemical potential plays 
the role of a "true" temperature, we have measured the 
function X(B) by quenching the system at different sub- 
critical values. Figure || shows that the slope of the FDT 
violating line, as well as the location of the breaking 
point, are independent of the quench value of at 
least in the time window explored in our numerical exper- 
iments and in the range of values of 1/fi used. Therefore 
as suggested in ref. jnj, the FDR appears locked in at 
the value the system would exhibit at the critical density, 
and should be only related with "universal" properties of 
the model, like for instance the exponent <f>. 

Conclusions. To summarize we have shown that three 
dimensional lattice-gas models defined by short range ki- 
netic constraint and trivial equilibrium Boltzmann-Gibbs 
measure exhibit a simple aging behavior in the integrated 
response function and violate fluctuation-dissipation the- 
orem in a manner similar to mean field models of struc- 
tural glasses. The fluctuation-dissipation ratio in the 
glassy phase appears to be a constant independent on 
the chemical potential of the reservoir; however, different 



results are expected if the boundary condition is changed 
or if creation and destruction of particles within the bulk 
is allowed. 

Given the non-holonomic nature of kinetic constraints 
and the trivial hamiltonian of the model, any statisti- 
cal mechanics approach based on the calculation of some 
partition function should be unsuitable to grasp the fea- 
tures of the glassy phase. In order to understand if aging 
is interrupted after a certain characteristic time [Q, it 
would be interesting to analytically investigate the long- 
time dynamical behavior of the model. A persistent ag- 
ing scenario would prove that also in finite-dimensional 
short-range models the glass transition may have a purely 
dynamic origin unrelated to any underlying thermody- 
namic singularity. 
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FIG. 2. Scaling behavior ol integrated response function 
X(t + t w ,t w ) vs. t/r w , where r w = t w + to with to = 10 3 
and t w = 10, 10 2 , 10 3 , 10 4 , 10 5 . The full line corresponds to 
X(t + t w ,t w ) ~ log(l + t/r w ). 
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1.5 2 

B(t+t„,t„) 

FIG. 3. Parametric plot of 2x(t + t w ,t w )/e vs. mean square 
displacement B(t + t w ,t w ) measured after a waiting time 
t w — 10 from a quench to the subcritical value 1/fi = 1/2.2. 
The density of initial configuration is 0.75. The two straight 
line have slope 1 (FDT line) and 0.79. Average over 200 sam- 
ples; error bars comparable to symbol size. 



FIG. 1. Integrated response function x(^ + ^>^) vs - time 
t after a quench to the subcritical values 1/fi — 1/2.2. The 
measurements were carried out after the perturbation was 
applied at waiting times t w — 10, 10 2 , 10 3 , 10 4 , 10°. The initial 
density is 0.75. Average over 20 samples. 
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B(t+t„,t„) 

FIG. 4. Parametric plot of 2x(t + t w ,t m )/e vs. mean square 
displacement B(t + t w ,t w ) measured after a waiting time 
t w = 10 5 from a quench to the subcritical value 1/fx = 1/2.05 
(stars), 1/2.1 (triangles), 1/2.2 (diamonds), 1/2.3 (squares). 
The density of initial configuration is 0.75. The two straight 
line have slope 1 and 0.79. Average over 200 samples; error 
bars comparable to symbols size. 
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